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Abstract 

The background in the region of the Q-value for neutrinoless double 
beta decay of 76 Ge has been investigated by different methods: Simu- 
lation with GEANT 4 of the HEIDELBERG-MOSCOW experiment, 
analysis of other Ge double beta experiments. Statistical features of 
the analysis at very low count rates are discussed. 



1 Introduction 



Recently first experimental evidence has been reported for neutrinoless dou- 
ble beta decay. Analysis of 55 kg y of data, taken by the HEIDELBERG- 
MOSCOW experiment in the GRAN SASSO over the years 1990 - 2000, has 
led [H El HI to a half-life 

T x /2 = (0.8 - 18.3) x 10 25 years (95% C.L.), (1) 

with best value of Ti/ 2 = 1.5 x 10 25 y, for the decay of the double beta emitter 
76 Ge 



76 Ge — > 76 Se + 2e~ (2) 

Assuming the decay amplitude to be dominated by exchange of a massive 
Majorana neutrino (see, e.g. [S]), this half-life results in a value of the effec- 
tive neutrino mass 

< m >= |£ [/2 TO .| = o.05 - 0.84 eV (95% C.L.), (3) 

with best value of 0.39 eV. Here a 50% uncertainty in the nuclear matrix 
elements has been taken into account (for details see j2J). 



This is for the first time that the absolute scale of the neutrino mass spectrum 
has been fixed, which cannot be achieved by neutrino oscillation experiments. 
This result restricts possible neutrino mass scenarios to degenerate or (still 
marginally allowed) inverse hierarchy |S1 El IH| • In the degenerate case it leads 
to a common neutrino mass eigenvalue of 

mx =0.05 -3.4 eV (95% C.L.). (4) 

This result is nicely consistent with later collected or analyzed experimental 
data, such as results from Large Scale Structure and CMB measurements 
[HJEI1E], or ultra-high energy cosmic rays The former yield an upper 
limit of Si m i=l-0 eV (corresponding in the degenerate case to a common 
mass eigenvalue mo < 0.33 eV). The Z-burst scenario for ultra-high energy 
cosmic rays requires 0.1 - 1.3 eV [T2*] . Tritium single beta decay cuts the 
upper range in eq.((U) down to 2.2 or 2.8 eV [T3]. 



There is further theoretical support for a neutrino mass in the range fixed by 
the HEIDELBERG-MOSCOW experiment. A model based on an A4 sym- 
metry of the neutrino mass matrix requires the neutrinos to be degenerate 
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Figure 1: The spectrum taken with 76 Ge detectors Nrs. 1,2,3,4,5 over the period 
August 1990 - May 2000 (54.98 kg y), in the energy range 2000 - 2100 keV. 
Simultaneous fit of the 214 Bi-lines and the two high-energy lines yield a probability 
for a line at 2039 keV of 91% C.L. 



and the common mass eigenvalue to be >0.2 eV |14j . 

Starting with the hypothesis that quark and lepton mixing are identical at 
or near the GUT scale, Mohapatra et. al. ^5] show that the large solar 
and atmospheric neutrino mixing angles can be understood purely as result 
of renormalization group evolution, if neutrino masses are quasi-degenerate 
(with same CP parity). The common Majorana neutrino mass then must be 
larger than 0.1 eV. 

In this paper we report additional support of the double beta result of 
[D El El E| , by further discussion of the structure of the experimental back- 
ground, by statistical considerations and by analysis of other double beta 
experiments investigating the decay of 76 Ge. 

Important points in the analysis of the measured spectrum are the structure 
of the background around Q^p (=2039.006(50) keV [IE!), an d the energy 
range of analysis around Qpp. 
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Figure 2: Result of the peak-searching procedure performed on the 
HEIDELBERG-MOSCOW spectrum using the Maximum Likelihood approach 
(left) and the Bayesian method (right). On the y axis the probability of having a 
line at the corresponding energy in the spectrum is shown. 




Figure 3: Result of the peak-search procedure performed for the UCBS/LBL 
spectrum |21j (left: Maximum Likelihood method, right: Bayes method). On the 
y axis the probability of having a line at the corresponding energy in the spectrum 
is shown. 



2 Structure of Background Around in 
Different Ge Double Beta Experiments 

Fig. 1 shows the spectrum measured in the range 2000 - 2100 keV in its orig- 
inal binning of 0.36 keV. By the peak search procedure developped (21 El on 
basis of the Bayes and Maximum Likelihood Methods, exploiting as impor- 
tant input parameters the experimental knowledge on the shape and width of 
lines in the spectrum, weak lines of 214 Bi have been identified at the energies 
of 2010.78, 2016.7, 2021.6 and 2052.94 keV P 12 El H • Fig. 2 shows the 
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Figure 4: Analysis of the spectrum measured by D. Caldwell et al. |21| . with 
the Maximum Likelihood Method, in the energy range 2000 - 2060 keV assuming 
lines at 2010.78, 2016.70, 2021.60, 2052.94, 2039.0 keV. No indication for a signal 
at 2039 kev is observed in this case. 



probability that there is a line of correct width and of Gaussian shape at a 
given energy, assuming all the rest of the spectrum as flat background (which 
is a highly conservative assumption), for the HEIDELBERG-MOSCOW ex- 
periment. 

The intensities of these lines have been shown to be consistent with other, 
strong Bi lines in the measured spectrum according to the branching ratios 
given in the Table of Isotopes ^Zj, and to Monte Carlo simulation of the 
experimental setup [3]. Note that the 2016 keV line, as an E0 transition, can 
be seen only by coincident summing of the two successive lines E = 1407.98 
keV and E = 609.316 keV. Its observation proves that the 238 U impurity 
from which it is originating, is located in the Cu cup of the detectors. Re- 
cent measurements of the spectrum of a 214 Bi source as function of distance 
source-detector confirm this interpretation [18 . 



Premature estimates of the Bi intensities given in Aalseth et.al,|hep-ex 



and Feruglio et.al., Nucl. Phys. B 637 (2002), 345, are incorrect, because this 
long-known spectroscopic effect of true coincident summing [2H] has not been 
taken into account, and also no simulation of the setup has been performed 
(for details see 0113) ■ 

These 214 Bi lines occur also in other investigations of double beta decay with 
Ge - and - even more important - also the additional structures in Fig. 2, 
which cannot be attributed at present, are seen in these other investigations. 

There are three other Ge experiments which have looked for double beta 
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decay of 76 Ge. First there is the experiment by Caldwell et al. [21], using 
natural Germanium detectors (7.8% abundance of 76 Ge, compared to 86% 
in the HEIDELBERG-MOSCOW experiment). This was the most sensitive 
natural Ge experiment. With their background a factor of 9 higher than in 
the HEIDELBERG-MOSCOW experiment and their measuring time of 22.6 
kg y ears, they had a statistics of the background by a factor of almost four 
larger than in the HEIDELBERG-MOSCOW experiment. This gives 
useful information on the composition of the background. 
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Figure 5: Result of the peak-search procedure performed for the ITEP/YePI 
spectrum |22j (left: Maximum Likelihood method, right: Bayes method). On the 
y axis the probability of having a line at the corresponding energy in the specrtum 
is shown. 
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Figure 6: Result of the peak-search procedure performed for the IGEX spectrum 
[SSlEl] using the ML approach (left) and the Bayesian statistics (right). On the y 
axis the probability of having a line at the corresponding energy in the spectrum 
is shown. 



Applying the same method of peak search as used in Fig. El yields indica- 
tions for peaks essentially at the same energies as in Fig. El (see Fig. This 
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shows that these peaks are not fluctuations. In particular it sees the 2010.78, 
2016.7, 2021.6 and 2052.94 keV 214 Bi lines, but a 1 s o the unattributed lines 
at higher energies. It finds, however, n o line at 2039 keV. This is consistent 
with the expectation from the rate found in the HEIDELBERG-MOSCOW 
experiment. About 16 observed events in the latter correspond to to 0.6 ex- 
pected events in the Caldwell experiment, because of the use of non-enriched 
material and the shorter measuring time. Fit of the Caldwell spectrum al- 
lowing for the 214 Bi lines and a 2039 keV line yields 0.4 events for the latter 
(see Fig. 0J. 

The first experiment using enriched (but not high-purity) Germanium 76 de- 
tectors was that of Kirpichnikov and coworkers [22] ■ These authors show only 
the energy range between 2020 and 2064 keV of their measured spectrum. 
The peak search procedure finds also here indications of lines around 2028 
keV and 2052 keV (see fig. but not any indication of a line at 2039 keV. 
This is consistent with the expectation, because for their low statistics of 
2.95 kg y they would expect here (according to HEIDELBERG-MOSCOW) 
0.9 counts. 

Another experiment (IGEX) used between 6 and 8.8 kg of enriched 76 Ge, 
but collected since beginning of the experiment in the early nineties till shut- 
down in 1999 only 8.8 kgyears of statistics [23 121- The authors of [21 121 
unfortunately show only the range 2020 to 2060 keV of their measured spec- 
trum in detail. Fig. [HI shows the result of our peak scanning of this range. 
Clear indications are seen for the Bi lines at 2021 and 2052 keV, but also of 
the unidentified structure around 2030 keV. Because of the conservative as- 
sumption on the background treatment in the scanning procedure (see above) 
there is no chance to see a signal at 2039 keV because of the 'hole' in the 
background of that spectrum (see Fig. 1 in [22]). With some good will one 
might see, however, an indication of 3 events here, consistent with the ex- 
pectation of the HEIDELBERG-MOSCOW experiment of 2.6 counts. 

It may be noted that all three experiments with enriched 76 Ge see a line 
around 2028 keV. It was already mentioned in [23] , that it is suspicious that 
this line which is seen also clearly in the pulse-shape analyzed spectrum of the 
HEIDELBERG-MOSCOW experiment (see [UEIE!) differs in energy from 
Qpp j us t by the K-shell X-ray energies of Ge (Se) of 9.2 (10.50) keV, or the 
K shell electron atomic binding energies of 11.10 (12.66) keV (see |H]), an d 
that this might point to a partly not understood mechanism of the double 
beta decay process. 
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Figure 7: Example of a random-generated spectrum with a Poisson distributed 
background with 4.0 events per channel and a Gaussian line centered in channel 
50 (line-width corresponds to a standard-deviation of a = 4.0 channels). The left 
picture shows a spectrum with a line-intensity of 10 events, the right spectrum a 
spectrum with a line-intensity of 100 events. The background is shown dark, the 
events of the line bright. 

3 Statistics at Low Count Rates - Peak Search, 
and Analysis Window 

At this point it may be useful to demonstrate the potential of the used peak 
search procedure. Fig. [7| shows a spectrum with Poisson-generated back- 
ground of 4 events per channel and a Gaussian line with width (standard 
deviation) of 4 channels centered at channel 50, with intensity of 10 (left) 
and 100 (right) events, respectively. Fig. El right shows the result of the anal- 
ysis of spectra of different line intensity with the Bayes method (here Bayes 
1-4 correspond to different choice of the prior distribution: (1) fi(r]) = 1 
(flat), (2) fi(rj) = 1/r), (3) fj,(r)) = l/y^rj, (4) Jeffrey's prior) and the Maxi- 
mum Likelihood Method. For each prior 1000 spectra have been generated 
with equal background and equal line intensity using random number gen- 
erators available at CERN [2"T| . The average values of the best values agree 
(see Fig. |HJ) very well with the known intensities also for very low count rates 
(as in Fig. left). 

In Fig. IHlwe show two simulations of a Gaussian line of 15 events, centered 
at channel 50, again with width (standard deviation) of 4 channels, on a 
Poisson-distributed background with 0.5 events/channel. The figure gives 
an indication of the possible degree of deviation of the energy of the peak 
maximum from the transition energy, on the level of statistics collected in ex- 
periments like the HEIDELBERG-MOSCOW experiment (here one channel 
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Figure 8: Results of analysis of random- number generated spectra, using Bayes 
and Maximum Likelihood method (the first one with different prior distributions). 
For every number of events in the simulated line, shown on the x-axis, 1000 random 
generated spectra were evaluated with the five given methods. The analysis on 
the left side was performed with an Poisson distributed background of 0.5 events 
per channel, the background for the spectra on the right side was 4.0 events per 
channel. Every vertical line shows the mean value of the calculated best values 
(thick points) with the la error area. The mean values are in good agreement 
with the expected values (black horizontal dashed lines). 

corresponds to 0.36 keV). This should be considered when comparing Figs. 

HE1E1E1 

The influence of the choice of the energy range of the analysis around Qpp 
has been thoroughly discussed in EJ- Since erroneous ideas about this 
point are still around, a few further comments may be given here. In Fig. El 
we show the analysis of a simulated spectrum consisting of a Gaussian line 
of width (standard deviation) of 4 channels and intensity of 10 counts on a 
Poisson-distributed background of 4 events per channel (see fig. [3 left) , as 
function of the width of the range of analysis. It is seen that a reliable result 
is obtained for a range of analysis of not smaller than 35 channels (i.e. ±18 
channels) - one channel corresponding to 0.36 keV in the HEIDELBERG- 
MOSCOW experiment. This is an important result, since it is of course 
important to keep the range of analysis as s m a 1 1 as possible, to avoid to 
include lines in the vicinity of the weak signal into the background. This 
unavoidably occurs when e.g. proceeding as suggested in Aalseth et. al., 
hep-ex/02020T8"l and Mod. Phys. Lett. A 17 (2002) 1475, Yu.G. Zdesenko 
et. al., Phys. Lett. B 546 (2002) 206. The arguments given in those papers 
are therefore incorrect. Also Kirpichnikov, who states j2H] that his analysis 
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Figure 9: Two spectra with a Poisson-distributed background and a Gaussian 
line with 15 events centered in channel 50 (with a width (standard-deviation) of 
4.0 channels) created with different random numbers. Shown is the result of the 
peak-scanning of the spectra. In the left picture the maximum of the probability 
corresponds well to the expected value (black line) whereas in the right picture a 
larger deviation is found. When a channel corresponds to 0.36 keV the deviation 
in the right picture is ~ 1.44 keV. 



finds a 2039 keV signal in the HEIDELBERG-MOSCOW spectrum on a 4 
sigma confidence level (as we also see it, when using the Feldman- Cousins 
method [26]), makes this mistake when analyzing the pulse-shape spectrum. 



4 Simulation of the Background with 
GEANT 4 

Finally the background around Qpp will be discussed from the side of simu- 
lation. A very careful new simulation of the different components of radioac- 
tive background in the HEIDELBERG-MOSCOW experiment has been per- 
formed recently by a new Monte Carlo program based on GEANT4 j2U EH] . 
This simulation uses a new event generator for simulation of radioactive de- 
cays basing on ENSDF-data and describes the decay of arbitrary radioactive 
isotopes including alpha, beta and gamma emission as well as conversion 
electrons and X-ray emission. Also included in the simulation is the influ- 
ence of neutrons in the energy range from thermal to high energies up to 
100 MeV on the measured spectrum. Elastic and inelastic reactions, and 
capture have been taken into account, and the corresponding production of 
radioactive isotopes in the setup. The neutron fluxes and energy distribu- 
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Figure 10: Result of an analysis as function of the evaluation width. The used 
spectrum consists of a Poisson distributed background with 4 events per channel, 
and a line of 10 events (see fig. left part). The dark area corresponds to a 68.3% 
confidence area with the dark line being the best value. Below an evaluation width 
of 35 channels the result becomes unreliable, above 35 channels the result is stable. 

tions were taken from published measurements performed in the Gran Sasso. 
Also simulated was the cosmic muon flux measured in the Gran Sasso, on 
the measured spectrum. To give a feeling for the quality of the simulation, 
Fig. ^2 shows the simulated and the measured spectra for a 228 Th source 
spectrum for 4 of our five detectors. The agreement is excellent. 

The simulation of the background of the experiment reproduces a 1 1 observed 
lines in the energy range between threshold (around 100 keV) and 2020 keV 
|2Sj. Fig. fT2l shows the simulated background in the range 2000-2100 keV 
with all k n o w n background components. The black solid line corresponds 
to the measured data in the period 20.11.1995 - 16.4.2002 (55.57 kg y). 

The background around Qpp is according to the simulations flat, the 
only expected lines come from 214 Bi (from the 238 U natural decay chain) at 
2010.89, 2016.7, 2021.6, 2052.94, 2085.1 and 2089.7 keV. Lines from cosmo- 
genically produced 56 Co (at 2034.76 keV and 2041.16 keV), half-life 77.3 days, 
are not expected since the first 200 days of measurement of each detector are 
not used in the data analysis. Also the potential contribution from decays 
of 77 Ge, 66 Ga, or 228 Ac, should not lead to signals visible in our measured 
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Figure 11: Comparison of the measured data (black line, November 1995 to April 
2002) and simulated spectrum (red line) for the detectors Nrs. 1,2,3 and 5 for a 
232 Th source spectrum. The agreement of simulation and measurement is excellent. 



spectrum near the signal at Qpp. For details we refer to [30] . 

The structures around 2028 keV, 2066 keV and 2075 keV seen - as also the 
214 Bi lines - in practically all Ge experiments (see above), cannot be identified 
at present. The 2028 keV line because of its strong occurence in the PSA 
spectrum in the HEIDELBERG-MOSCOW experiment, may play a special 
role here (see above). 



5 Conclusion 

Concluding, additional support has been given for the evidence of a signal 
for neutrinoless double beta decay, by showing consistency of the result - for 
the signal, a n d for the background - with other double beta decay exper- 
iments using non-enriched or enriched Germanium detectors. In particular 
it has been shown that the lines seen in the vicinity of the signal (including 
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Figure 12: Simulated background of the HEIDELBERG-MOSCOW experiment 
in the energt range from 2000 to 2100 keV with all known background components. 
The black trough-drawn line corresponds to the measured data from 20.11.1995 to 
16.4.2002 (55.57 kg y). 



those which at present cannot be attributed) are seen also in the other ex- 
periments. This is important for the correct treatment of the background. 
Furthermore, the sensitivity of the peak identification procedures has been 
demonstrated by extensive statistical simulations. It has been further shown 
by new extensive simulations of the expected background by GEANT4, that 
the background around Qpp should be flat, and that no known gamma line 
is expected at the energy of Qpp. The 2039 keV signal is seen only in 
the HEIDELBERG-MOSCOW experiment, which has a factor of 10, and 
much more, statistics than all other double beta experiments. 
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Abstract 

The background in the region of the Q-value for neutrinoless double 
beta decay of 76 Ge has been investigated by different methods: Simu- 
lation with GEANT 4 of the HEIDELBERG-MOSCOW experiment, 
analysis of other Ge double beta experiments. Statistical features of 
the analysis at very low count rates are discussed. 



1 Introduction 



Recently first experimental evidence has been reported for neutrinoless dou- 
ble beta decay. Analysis of 55 kg y of data, taken by the HEIDELBERG- 
MOSCOW experiment in the GRAN SASSO over the years 1990 - 2000, has 
led [H El HI to a half-life 

T x /2 = (0.8 - 18.3) x 10 25 years (95% C.L.), (1) 

with best value of Ti/ 2 = 1.5 x 10 25 y, for the decay of the double beta emitter 
76 Ge 



76 Ge — > 76 Se + 2e~ (2) 

Assuming the decay amplitude to be dominated by exchange of a massive 
Majorana neutrino (see, e.g. [S]), this half-life results in a value of the effec- 
tive neutrino mass 

< m >= |£ [/2 TO .| = o.05 - 0.84 eV (95% C.L.), (3) 

with best value of 0.39 eV. Here a 50% uncertainty in the nuclear matrix 
elements has been taken into account (for details see j2J). 



This is for the first time that the absolute scale of the neutrino mass spectrum 
has been fixed, which cannot be achieved by neutrino oscillation experiments. 
This result restricts possible neutrino mass scenarios to degenerate or (still 
marginally allowed) inverse hierarchy |S1 El IH| • In the degenerate case it leads 
to a common neutrino mass eigenvalue of 

mx =0.05 -3.4 eV (95% C.L.). (4) 

This result is nicely consistent with later collected or analyzed experimental 
data, such as results from Large Scale Structure and CMB measurements 
[HJEI1E], or ultra-high energy cosmic rays The former yield an upper 
limit of Si m i=l-0 eV (corresponding in the degenerate case to a common 
mass eigenvalue mo < 0.33 eV). The Z-burst scenario for ultra-high energy 
cosmic rays requires 0.1 - 1.3 eV [T2*] . Tritium single beta decay cuts the 
upper range in eq.((U) down to 2.2 or 2.8 eV [T3]. 



There is further theoretical support for a neutrino mass in the range fixed by 
the HEIDELBERG-MOSCOW experiment. A model based on an A4 sym- 
metry of the neutrino mass matrix requires the neutrinos to be degenerate 
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Figure 1: The spectrum taken with 76 Ge detectors Nrs. 1,2,3,4,5 over the period 
August 1990 - May 2000 (54.98 kg y), in the energy range 2000 - 2100 keV. 
Simultaneous fit of the 214 Bi-lines and the two high-energy lines yield a probability 
for a line at 2039 keV of 91% C.L. 



and the common mass eigenvalue to be >0.2 eV |14j . 

Starting with the hypothesis that quark and lepton mixing are identical at 
or near the GUT scale, Mohapatra et. al. ^5] show that the large solar 
and atmospheric neutrino mixing angles can be understood purely as result 
of renormalization group evolution, if neutrino masses are quasi-degenerate 
(with same CP parity). The common Majorana neutrino mass then must be 
larger than 0.1 eV. 

In this paper we report additional support of the double beta result of 
[D El El E| , by further discussion of the structure of the experimental back- 
ground, by statistical considerations and by analysis of other double beta 
experiments investigating the decay of 76 Ge. 

Important points in the analysis of the measured spectrum are the structure 
of the background around Q^p (=2039.006(50) keV [IE!), an d the energy 
range of analysis around Qpp. 
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Figure 2: Result of the peak-searching procedure performed on the 
HEIDELBERG-MOSCOW spectrum using the Maximum Likelihood approach 
(left) and the Bayesian method (right). On the y axis the probability of having a 
line at the corresponding energy in the spectrum is shown. 




Figure 3: Result of the peak-search procedure performed for the UCBS/LBL 
spectrum |21j (left: Maximum Likelihood method, right: Bayes method). On the 
y axis the probability of having a line at the corresponding energy in the spectrum 
is shown. 



2 Structure of Background Around in 
Different Ge Double Beta Experiments 

Fig. 1 shows the spectrum measured in the range 2000 - 2100 keV in its orig- 
inal binning of 0.36 keV. By the peak search procedure developped (21 El on 
basis of the Bayes and Maximum Likelihood Methods, exploiting as impor- 
tant input parameters the experimental knowledge on the shape and width of 
lines in the spectrum, weak lines of 214 Bi have been identified at the energies 
of 2010.78, 2016.7, 2021.6 and 2052.94 keV P 12 El H • Fig. 2 shows the 
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Figure 4: Analysis of the spectrum measured by D. Caldwell et al. |21| . with 
the Maximum Likelihood Method, in the energy range 2000 - 2060 keV assuming 
lines at 2010.78, 2016.70, 2021.60, 2052.94, 2039.0 keV. No indication for a signal 
at 2039 kev is observed in this case. 



probability that there is a line of correct width and of Gaussian shape at a 
given energy, assuming all the rest of the spectrum as flat background (which 
is a highly conservative assumption), for the HEIDELBERG-MOSCOW ex- 
periment. 

The intensities of these lines have been shown to be consistent with other, 
strong Bi lines in the measured spectrum according to the branching ratios 
given in the Table of Isotopes ^Zj, and to Monte Carlo simulation of the 
experimental setup [3]. Note that the 2016 keV line, as an E0 transition, can 
be seen only by coincident summing of the two successive lines E = 1407.98 
keV and E = 609.316 keV. Its observation proves that the 238 U impurity 
from which it is originating, is located in the Cu cup of the detectors. Re- 
cent measurements of the spectrum of a 214 Bi source as function of distance 
source-detector confirm this interpretation [18 . 



Premature estimates of the Bi intensities given in Aalseth et.al,|hep-ex 



and Feruglio et.al., Nucl. Phys. B 637 (2002), 345, are incorrect, because this 
long-known spectroscopic effect of true coincident summing [2H] has not been 
taken into account, and also no simulation of the setup has been performed 
(for details see 0113) ■ 

These 214 Bi lines occur also in other investigations of double beta decay with 
Ge - and - even more important - also the additional structures in Fig. 2, 
which cannot be attributed at present, are seen in these other investigations. 

There are three other Ge experiments which have looked for double beta 
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decay of 76 Ge. First there is the experiment by Caldwell et al. [21], using 
natural Germanium detectors (7.8% abundance of 76 Ge, compared to 86% 
in the HEIDELBERG-MOSCOW experiment). This was the most sensitive 
natural Ge experiment. With their background a factor of 9 higher than in 
the HEIDELBERG-MOSCOW experiment and their measuring time of 22.6 
kg y ears, they had a statistics of the background by a factor of almost four 
larger than in the HEIDELBERG-MOSCOW experiment. This gives 
useful information on the composition of the background. 
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Figure 5: Result of the peak-search procedure performed for the ITEP/YePI 
spectrum |22j (left: Maximum Likelihood method, right: Bayes method). On the 
y axis the probability of having a line at the corresponding energy in the specrtum 
is shown. 
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Figure 6: Result of the peak-search procedure performed for the IGEX spectrum 
[SSlEl] using the ML approach (left) and the Bayesian statistics (right). On the y 
axis the probability of having a line at the corresponding energy in the spectrum 
is shown. 



Applying the same method of peak search as used in Fig. El yields indica- 
tions for peaks essentially at the same energies as in Fig. El (see Fig. This 
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shows that these peaks are not fluctuations. In particular it sees the 2010.78, 
2016.7, 2021.6 and 2052.94 keV 214 Bi lines, but a 1 s o the unattributed lines 
at higher energies. It finds, however, n o line at 2039 keV. This is consistent 
with the expectation from the rate found in the HEIDELBERG-MOSCOW 
experiment. About 16 observed events in the latter correspond to to 0.6 ex- 
pected events in the Caldwell experiment, because of the use of non-enriched 
material and the shorter measuring time. Fit of the Caldwell spectrum al- 
lowing for the 214 Bi lines and a 2039 keV line yields 0.4 events for the latter 
(see Fig. 0J. 

The first experiment using enriched (but not high-purity) Germanium 76 de- 
tectors was that of Kirpichnikov and coworkers [22] ■ These authors show only 
the energy range between 2020 and 2064 keV of their measured spectrum. 
The peak search procedure finds also here indications of lines around 2028 
keV and 2052 keV (see fig. but not any indication of a line at 2039 keV. 
This is consistent with the expectation, because for their low statistics of 
2.95 kg y they would expect here (according to HEIDELBERG-MOSCOW) 
0.9 counts. 

Another experiment (IGEX) used between 6 and 8.8 kg of enriched 76 Ge, 
but collected since beginning of the experiment in the early nineties till shut- 
down in 1999 only 8.8 kgyears of statistics [23 121- The authors of [21 121 
unfortunately show only the range 2020 to 2060 keV of their measured spec- 
trum in detail. Fig. [HI shows the result of our peak scanning of this range. 
Clear indications are seen for the Bi lines at 2021 and 2052 keV, but also of 
the unidentified structure around 2030 keV. Because of the conservative as- 
sumption on the background treatment in the scanning procedure (see above) 
there is no chance to see a signal at 2039 keV because of the 'hole' in the 
background of that spectrum (see Fig. 1 in [22]). With some good will one 
might see, however, an indication of 3 events here, consistent with the ex- 
pectation of the HEIDELBERG-MOSCOW experiment of 2.6 counts. 

It may be noted that all three experiments with enriched 76 Ge see a line 
around 2028 keV. It was already mentioned in [23] , that it is suspicious that 
this line which is seen also clearly in the pulse-shape analyzed spectrum of the 
HEIDELBERG-MOSCOW experiment (see [UEIE!) differs in energy from 
Qpp j us t by the K-shell X-ray energies of Ge (Se) of 9.2 (10.50) keV, or the 
K shell electron atomic binding energies of 11.10 (12.66) keV (see |H]), an d 
that this might point to a partly not understood mechanism of the double 
beta decay process. 
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Figure 7: Example of a random-generated spectrum with a Poisson distributed 
background with 4.0 events per channel and a Gaussian line centered in channel 
50 (line-width corresponds to a standard-deviation of a = 4.0 channels). The left 
picture shows a spectrum with a line-intensity of 10 events, the right spectrum a 
spectrum with a line-intensity of 100 events. The background is shown dark, the 
events of the line bright. 

3 Statistics at Low Count Rates - Peak Search, 
and Analysis Window 

At this point it may be useful to demonstrate the potential of the used peak 
search procedure. Fig. [7| shows a spectrum with Poisson-generated back- 
ground of 4 events per channel and a Gaussian line with width (standard 
deviation) of 4 channels centered at channel 50, with intensity of 10 (left) 
and 100 (right) events, respectively. Fig. El right shows the result of the anal- 
ysis of spectra of different line intensity with the Bayes method (here Bayes 
1-4 correspond to different choice of the prior distribution: (1) fi(r]) = 1 
(flat), (2) fi(rj) = 1/r), (3) fj,(r)) = l/y^rj, (4) Jeffrey's prior) and the Maxi- 
mum Likelihood Method. For each prior 1000 spectra have been generated 
with equal background and equal line intensity using random number gen- 
erators available at CERN [2"T| . The average values of the best values agree 
(see Fig. |HJ) very well with the known intensities also for very low count rates 
(as in Fig. left). 

In Fig. IHlwe show two simulations of a Gaussian line of 15 events, centered 
at channel 50, again with width (standard deviation) of 4 channels, on a 
Poisson-distributed background with 0.5 events/channel. The figure gives 
an indication of the possible degree of deviation of the energy of the peak 
maximum from the transition energy, on the level of statistics collected in ex- 
periments like the HEIDELBERG-MOSCOW experiment (here one channel 
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Figure 8: Results of analysis of random- number generated spectra, using Bayes 
and Maximum Likelihood method (the first one with different prior distributions). 
For every number of events in the simulated line, shown on the x-axis, 1000 random 
generated spectra were evaluated with the five given methods. The analysis on 
the left side was performed with an Poisson distributed background of 0.5 events 
per channel, the background for the spectra on the right side was 4.0 events per 
channel. Every vertical line shows the mean value of the calculated best values 
(thick points) with the la error area. The mean values are in good agreement 
with the expected values (black horizontal dashed lines). 

corresponds to 0.36 keV). This should be considered when comparing Figs. 

HE1E1E1 

The influence of the choice of the energy range of the analysis around Qpp 
has been thoroughly discussed in EJ- Since erroneous ideas about this 
point are still around, a few further comments may be given here. In Fig. El 
we show the analysis of a simulated spectrum consisting of a Gaussian line 
of width (standard deviation) of 4 channels and intensity of 10 counts on a 
Poisson-distributed background of 4 events per channel (see fig. [3 left) , as 
function of the width of the range of analysis. It is seen that a reliable result 
is obtained for a range of analysis of not smaller than 35 channels (i.e. ±18 
channels) - one channel corresponding to 0.36 keV in the HEIDELBERG- 
MOSCOW experiment. This is an important result, since it is of course 
important to keep the range of analysis as s m a 1 1 as possible, to avoid to 
include lines in the vicinity of the weak signal into the background. This 
unavoidably occurs when e.g. proceeding as suggested in Aalseth et. al., 
hep-ex/02020T8"l and Mod. Phys. Lett. A 17 (2002) 1475, Yu.G. Zdesenko 
et. al., Phys. Lett. B 546 (2002) 206. The arguments given in those papers 
are therefore incorrect. Also Kirpichnikov, who states j2H] that his analysis 
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Figure 9: Two spectra with a Poisson-distributed background and a Gaussian 
line with 15 events centered in channel 50 (with a width (standard-deviation) of 
4.0 channels) created with different random numbers. Shown is the result of the 
peak-scanning of the spectra. In the left picture the maximum of the probability 
corresponds well to the expected value (black line) whereas in the right picture a 
larger deviation is found. When a channel corresponds to 0.36 keV the deviation 
in the right picture is ~ 1.44 keV. 



finds a 2039 keV signal in the HEIDELBERG-MOSCOW spectrum on a 4 
sigma confidence level (as we also see it, when using the Feldman- Cousins 
method [26]), makes this mistake when analyzing the pulse-shape spectrum. 



4 Simulation of the Background with 
GEANT 4 

Finally the background around Qpp will be discussed from the side of simu- 
lation. A very careful new simulation of the different components of radioac- 
tive background in the HEIDELBERG-MOSCOW experiment has been per- 
formed recently by a new Monte Carlo program based on GEANT4 j2U EH] . 
This simulation uses a new event generator for simulation of radioactive de- 
cays basing on ENSDF-data and describes the decay of arbitrary radioactive 
isotopes including alpha, beta and gamma emission as well as conversion 
electrons and X-ray emission. Also included in the simulation is the influ- 
ence of neutrons in the energy range from thermal to high energies up to 
100 MeV on the measured spectrum. Elastic and inelastic reactions, and 
capture have been taken into account, and the corresponding production of 
radioactive isotopes in the setup. The neutron fluxes and energy distribu- 
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Figure 10: Result of an analysis as function of the evaluation width. The used 
spectrum consists of a Poisson distributed background with 4 events per channel, 
and a line of 10 events (see fig. left part). The dark area corresponds to a 68.3% 
confidence area with the dark line being the best value. Below an evaluation width 
of 35 channels the result becomes unreliable, above 35 channels the result is stable. 

tions were taken from published measurements performed in the Gran Sasso. 
Also simulated was the cosmic muon flux measured in the Gran Sasso, on 
the measured spectrum. To give a feeling for the quality of the simulation, 
Fig. ^2 shows the simulated and the measured spectra for a 228 Th source 
spectrum for 4 of our five detectors. The agreement is excellent. 

The simulation of the background of the experiment reproduces a 1 1 observed 
lines in the energy range between threshold (around 100 keV) and 2020 keV 
|2Sj. Fig. fT2l shows the simulated background in the range 2000-2100 keV 
with all k n o w n background components. The black solid line corresponds 
to the measured data in the period 20.11.1995 - 16.4.2002 (55.57 kg y). 

The background around Qpp is according to the simulations flat, the 
only expected lines come from 214 Bi (from the 238 U natural decay chain) at 
2010.89, 2016.7, 2021.6, 2052.94, 2085.1 and 2089.7 keV. Lines from cosmo- 
genically produced 56 Co (at 2034.76 keV and 2041.16 keV), half-life 77.3 days, 
are not expected since the first 200 days of measurement of each detector are 
not used in the data analysis. Also the potential contribution from decays 
of 77 Ge, 66 Ga, or 228 Ac, should not lead to signals visible in our measured 
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Figure 11: Comparison of the measured data (black line, November 1995 to April 
2002) and simulated spectrum (red line) for the detectors Nrs. 1,2,3 and 5 for a 
232 Th source spectrum. The agreement of simulation and measurement is excellent. 



spectrum near the signal at Qpp. For details we refer to [30] . 

The structures around 2028 keV, 2066 keV and 2075 keV seen - as also the 
214 Bi lines - in practically all Ge experiments (see above), cannot be identified 
at present. The 2028 keV line because of its strong occurence in the PSA 
spectrum in the HEIDELBERG-MOSCOW experiment, may play a special 
role here (see above). 



5 Conclusion 

Concluding, additional support has been given for the evidence of a signal 
for neutrinoless double beta decay, by showing consistency of the result - for 
the signal, a n d for the background - with other double beta decay exper- 
iments using non-enriched or enriched Germanium detectors. In particular 
it has been shown that the lines seen in the vicinity of the signal (including 
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Figure 12: Simulated background of the HEIDELBERG-MOSCOW experiment 
in the energt range from 2000 to 2100 keV with all known background components. 
The black trough-drawn line corresponds to the measured data from 20.11.1995 to 
16.4.2002 (55.57 kg y). 



those which at present cannot be attributed) are seen also in the other ex- 
periments. This is important for the correct treatment of the background. 
Furthermore, the sensitivity of the peak identification procedures has been 
demonstrated by extensive statistical simulations. It has been further shown 
by new extensive simulations of the expected background by GEANT4, that 
the background around Qpp should be flat, and that no known gamma line 
is expected at the energy of Qpp. The 2039 keV signal is seen only in 
the HEIDELBERG-MOSCOW experiment, which has a factor of 10, and 
much more, statistics than all other double beta experiments. 
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